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Abstract. The low-energy properties of two-dimensional ensembles of dipole-coupled magnetic 
nanoparticles are studied as function of structural disorder and particle coverage. Already small 
deviations from a square particle arrangement lift the degeneracies of the microvortex magnetic 
configuration, and result in a strongly noncollinear magnetic order of the particle ensemble. The 
energy distribution of metastable states is determined. For a low degree of disorder a strongly 
asymmetric shape with a pronounced peak of the ground state energy results. In contrast, for a 
strong disorder a Gaussian-like distribution is obtained. The average dipole coupling energy E^ip 
decreases with increasing structural disorder. The role of vacancies has been studied for a square 
particle array by determining the angular distribution of the preferred microvortex angle as function 
of the vacancy concentration. Indications for a preferred angular direction along the axial as well 
as along the diagonal directions of the square array are revealed. A corresponding investigation for 
disturbed square arrays results in a different angular distribution. The effect of dipole-quadrupole 
corrections resulting from the finite size of the particles is quantified. 
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1. Introduction 

Interacting magnetic nanostructured materials are currently the subject of intense research activity, 
driven by their fundamental interest and technological perspectives [1] . Numerous experimental and 
theoretical studies have been performed for various two-dimensional (2D) [2] and three-dimensional 
(3D) [3, 4, 5, 6, 7, 8] arrangements of nanometer-size magnetic particles having different degrees of 
structural and magnetic disorder. The magnetic behaviour of these systems is determined by single- 
particle properties (e.g., particle moments, lattice and shape anisotropies, etc.), by the composition 
and morphology of the nanostructure, and in particular by the nature of the dominant interparticle 
interactions. The latter comprises especially the magnetic dipole coupling, which will be adressed in 
the present study. Other interparticle interactions are, for example, the Rudermann-Kittel-Kasuya- 
Yosida (RKKY) indirect exchange mediated by the conduction electrons of a metallic substrate, or 
the short-range direct exchange in case when the particles are in contact. The relative importance of 
single-particle versus interparticle contributions can be tuned experimentally at least to some extent 
by changing sample characteristics such as the particle-size distribution or the average interparticle 
distance. For low particle coverages the interactions can be treated as a perturbation to the single- 
particle properties. However, for dense particle ensembles the interparticle interactions become 
increasingly important and eventually dominate. In this interesting case the single-particle approach 
is no longer applicable and an explicit treatment of the interactions is unavoidable [3, 4, 5]. 

A fundamental question in this context is to identify and understand the collectively ordered 
magnetic states which are induced by the interactions in such particle ensembles. A few basic 
properties of strictly periodic dipole-coupled systems are summarized. The ground state of a 
square lattice of equal-sized (monodispersed) particles is the so-called microvortex (MV) magnetic 
arrangement [9, 10]. This magnetic structure is characterized by the microvortex angle mv , where 
the angles of the particle magnetizations of a plaquette of four neighbouring particles are given 
by 0i = </W) 02 = ~ 0mv, 03 = 180° + mv , and 04 = 180° — mv . In particular, multiples of 
0mv = 90° represent columnar states, consisting of ferromagnetic rows or columns with alternating 
signs of magnetizations. Evidently, the MV state has a vanishing net magnetization. The parallel 
or ferromagnetic (FM) state has a larger energy and is actually an unstable solution. The ordering 
of dipole-coupled spins is very sensitive to the lattice structure. For example, the honeycomb lattice 
has a ground state with a vanishing net magnetization which is similar to the MV state of the square 
lattice [10]. In contrast, for the hexagonal lattice the ground state is ferromagnetic [11, 12]. 

Despite the fact that the dipole interaction is not rotationally invariant, the ground state of 
these dipole-coupled periodic particle arrays are continuously degenerate with respect to a rotation 
of the MV angle mv or of the FM angle fm , respectively. This holds for classical spins at 
T = 0. Thermal fluctuations, quantum fluctuations, or a structural disorder immediately lift these 
accidental ground state degeneracies of the periodic structures. The energy lowering associated to 
the symmetry breaking stabilizes the system in some particular magnetic order. Therefore, the 
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square and honeycomb lattices are said to show a so-called 'order-by-disorder' effect [10, 13]. This 
implies that internal (nontrivial) degeneracies (e.g., the relative directions of magnetic sublattices) 
are lifted by the presence of disorder, and that particular directions of the sublattice magnetizations 
are preferred. Moreover, for a dipole-coupled square lattice the fluctuations induce a fourfold in- 
plane magnetic anisotropy. As shown by Prakash and Henley [10], the preferred in-plane magnetic 
orientations are the axial directions (columnar states) for thermal and quantum fluctuations, and the 
diagonal directions for a structural disorder induced by a small amount of vacancies. In addition, 
Monte Carlo calculations and interacting spin wave theory indicate that a magnetic ordering at finite 
temperatures exists for a dipole-coupled square spin lattice, since the magnetic excitations are not 
continuously degenerate [14]. Thus, the Mermin- Wagner theorem, which excludes an ordered state 
for a continuously degenerate 2D system at finite temperatures [15], is not effective in this case. 
Notice that already the long-range character of the dipole interaction induces a collectively ordered 
magnetic state in a square lattice [16]. 

Structural disorder, which is usually present in real magnetic nanostructures, results for example 
from the size and shape dispersion, positional disorder, or random anisotropy axes. Due to the 
nonuniform and competing nature of the magnetic couplings the magnetic ordering in disordered 
particle ensembles is similar to the one of a spin-glass [17]. Thus, many different metastable 
states exist, which are characterized by strong magnetic noncollinearities. An intriguing question 
in this context is whether the dipole interaction in disordered planar particle ensembles results in 
a collectively ordered magnetic state. Several experiments on interacting, high-density ferrofluid 
systems indicate the onset of a collectively ordered state below a characteristic, concentration- 
dependent temperature. For example, at this temperature a 'critical slowing down' of the magnetic 
relaxation is observed [3, 4]. Furthermore, recent measurements on Co islands on Cu(001) exhibit a 
magnetic hysteresis and remanence in the temperature range up to 150 K also for coverages below 
the magnetic percolation threshold [18]. Due to the small size of the Co islands these findings could 
not explained simply by single-particle blocking effects. Note that the experimental determination 
of the ordering temperature is difficult, since the relaxation times are often very long. Hence, it is 
of considerable interest to analyze how the magnetic order depends on the sample parameters which 
can be controlled in experiment. 

The purpose of this paper is the theoretical study of the low-energy properties of dipole-coupled 
magnetic particle ensembles forming inhomogeneous planar arrangements with various degrees of 
disorder ranging from a quasi-periodic square lattice to a random array. Numerical simulations are 
performed in order to achieve a detailed microscopic description within the model. The low symmetry 
of the system and the complicated nature of the interaction seem to preclude simple analytical 
approaches. We focus on the strongly interacting case, neglecting single-particle anisotropies. Special 
attention is paid to the role of disorder and noncollinear arrangements of the particle magnetic 
moments. In particular, we determine the energy distribution of metastable states and the average 
magnetic dipole energy for different coverages and types of disorder. Furthermore, we investigate 
global and local order parameters as function of disorder, and compare our results with previous 
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calculations [10]. All calculations are performed at T = 0. 

The rest of the paper is organized as follows. The theoretical methods are outlined in section 2. 
Representative results for the magnetic arrangements, as well as for various magnetic properties are 
presented in section 3. A conclusion is given in section 4. Finally, in the Appendix we report details 
for the extension of the dipole summation beyond the point-dipole approximation which takes into 
account effects resulting from the finite particle size. 



2. Theory 



We consider a 2D rectangular unit cell in the xy- plane with n = n x x n y non-overlapping, disk- 
shaped magnetic particles. Due to the strong direct exchange interaction and the small size of the 
particles under consideration, each particle i can be viewed as a single magnetic domain (Stoner- 
Wohlfarth particle) [19]. Thus, a particle containing Ni atoms carries a giant spin = Ni/i at , 
where /J, at is the atomic magnetic moment. For simplicity, we restrict the particle magnetizations 
to be confined to the xy plane, the planar rotator Mi is then characterized by the in-plane angle 
(f)f Mi = (Mf,Mf,M[) = Mj(cos 0j, sin 0j, 0). In this study no size dispersion is considered, 
i.e., Ni = N. For disk-shaped particles the particle radius r is given by r /a = y/N, where 
a is the interatomic distance. Unless otherwise stated, the present results refer to a unit cell 
containing n = 100 particles with N = 1000 atoms each. The size of the unit cell is given by 
L x x L y = (n x Rq) x (n y R ), with Rq the average interparticle distance. For a planar array of 
circular particles the overall surface coverage is C = TT(r /R ) 2 . Four different types of the lateral 
particle arrangement have been considered: (i) a periodic square array, i.e., the particle centers are 
located on the sites of a square lattice with lattice constant Rq, (ii) a disturbed (quasi-periodic) array 
for which the particle centers deviate randomly from the square array, using a Gaussian distribution 
P(R) with positional standard deviation ctr, (iii) a diluted square particle lattice containing a number 
of vacancies n vac with concentration C vac , and (iv) a fully random distribution of non-overlapping 
particles within the unit cell. Periodic boundary conditions are introduced in order to describe an 
infinitely extended planar particle ensemble. By considering a single unit cell a finite system can be 
realized as well. In this case the effect of boundaries will strongly dominate the resulting magnetic 
arrangement. 

For such particle ensembles we consider the dipole-dipole interaction between the magnetic 
moments M i: which can also be expressed in terms of a dipole field Bf ip acting on Mi due to all 
other particle magnetic moments: 

E ^ = yE[ M ' M > r ^ ~ 3fo M i) Mi) ^] =-yE Mi ^ P » (!) 

i,3 i 

where r^ = jr^j = | — r j- 1 is the distance between the centers of particles i and j, and fi the vacuum 
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permeability. The infinite range of the dipole interaction is taken into account by applying an Ewald- 
type summation over all periodically arranged unit cells of the extended thin film [20]. In addition 
to the usual point-dipole sum we consider the leading correction resulting from the finite particle size 
(dipole-quadrupole interaction), which is outlined in detail in the Appendix. This correction becomes 
comparable to the point-dipole sum for large particle coverages or small interparticle distances. The 
energy unit of the dipole coupling is given by E® ip = /j% t /a,Q, with ji at in units of the Bohr magneton 
li B . In this study we assume values appropriate to Fe (fj, at = 2.2 /i B , a = 2.5 A), which yields 
= 0.19 K. The effects of single-particle anisotropies resulting from the spin-orbit interaction, the 
dipole interaction among atomic magnetic moments within each particle (shape anisotropy), external 
magnetic fields, and finite temperatures are beyond the scope of the present study. 

Starting from an arbitrary initial configuration {0J mtia1 } of the magnetic directions, the total 
magnetic energy E dip of the system is relaxed to the nearest local minimum, which often corresponds 
to a metastable state, by varying all in-plane angles fa of the particles using a conjugated gradient 
method [21]. For example, the experimental situation of a remanent state after removal of an 
external magnetic field is simulated by choosing a fully aligned initial arrangement along a certain 
direction. We emphasize that the applied procedure is not intended to search preferently for the 
global energy minimum or ground state, which is the equilibrium state at T = 0. Rather, at first 
we determine the energy distribution of the local minima for different degrees of structural disorder 
introduced in the planar particle array. A number of randomly chosen intital configurations {$ mtia1 } 
is created and relaxed to a nearby local minimum. This metastable state is characterized by its 
energy E d i P and by its set of angles {<fii}. A twofold (uniaxial) symmetry is always present due 
to time inversion symmetry, i.e., the energy of a state does not change under the transformation 
0i - ► fa + 180° performed simultaneously for all particles. These mirrored states are considered 
to be equivalent. The numbers of trials yielding relaxed states with energies falling into given 
energy intervals are monitored, and the corresponding energy histogram is determined. Since the 
minimization procedure relaxes typically to the local energy minimum that is closest to the initial 
state, this frequency provides a measure of the catchment area, i.e., the area of the 'energy valley' 
belonging to that state in the high-dimensional configuration space. 

In order to account for the large number of local energy minima occurring in particle 
arrangements with structural inhomogeneities we determine the average dipole energy E d i P resulting 
from many different trials for the same particle array. In addition, we average -Edip over a number of 
different realizations of the unit cell, using the same global variables which characterize the particle 
ensemble (coverage, standard deviation a R , etc.). For comparison, E dip is also calculated for the same 
spatial setup by assuming a ferromagnetic state (i.e., fa = fa m ) and averaging over the FM angle 
fa m . In case of a disturbed (quasi-periodic) square lattice, i.e., if the particle numbering allows for 
the identification of a square plaquette of four neighbouring particles, we also determine the dipole 
energy for the microvortex magnetic arrangement, averaging over the MV angle <f) mv . Note that for a 
nonuniform particle ensemble these FM and MV magnetic arrangements do not correspond in general 
to local energy minima. For a random set of angles {0j} one obtains -Edip — which constitutes the 
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natural energy reference. 

The deviations of the metastable low-energy magnetic arrangements from the MV state are 
quantified by the global and local microvortex order parameters M|J° bal (o"R) and M^^ctr) given by 

1 r / \ 2 / \ 2 i 1 / 2 

M£*(vr) = -[(E^ 1 )** 008 ^) + (E(- 1 )" sin ^) ] > ( 2 ) 

i i 

M =^[(E(-ir a cos0 i ) 2 +(E(-ir^0,) 2 ] 1/2 , (3) 



s^local ' ' 

mv 



where % y and i x denote the numbering of the rows and columns corresponding to particle %. Clearly, 
a reference square lattice is prerequisitive. Hence, we restrict ourselves to a quasi-periodic array 
or to a square array with vacancies for which i y and i x can be uniquely defined. These two order 
parameters differ by the sum % running either over all n particles of the unit cell or over the four 
neighbouring particles within a square plaquette. They have a simple geometrical interpretation as 
the projection of the magnetic configuration on the two linear independent columnar states having 
0mv = 0° and mv = 90°, regardless of all possible rotations of the MV states within this plane. 
M m 1 ° bal (o"R) measures the MV order in the nanostructure, whereas M^^or) measures the short 
range ordering. Both quantities are averaged over an appropriate number of initial configurations for 
the same particle arrangement and over different realizations of the unit cell. In addition M^ al (<7jj) 
is averaged over all four-particle plaquettes within the unit cell. 

Finally, another quantity characterizing the magnetic properties of nanostructured particle 
arrangements is the distribution of the microvortex angles mv . Prakash and Henley observed that 
the preferred MV angles for a square particle array with a small number of randomly distributed 
vacancies are the diagonal directions (0 mv = 45°, 135°, etc.) [10]. It is therefore interesting to 
investigate how the angular distribution of the microvortex angles depends on the type and degree of 
disorder. For this purpose different realizations of the nanostructure are created. For each particle 
setup the MV state with angle mv yielding the lowest energy is determined, and the corresponding 
frequencies for all angles are monitored. We consider a disturbed square particle array with standard 
deviation ctr, and a diluted square particle lattice with a concentration C vac of vacancies. It should 
be noted that for a small number of vacancies one obtains an appreciable dependence on the shape of 
the unit cell (aspect ratio of the rectangle), since these vacancies are located effectively on a square 
or rectangular lattice due to the periodic boundary conditions. 



3. Results 



First, figure 1 illustrates some representative low-energy magnetic arrangements of disordered particle 
arrays. With increasing disorder, characterized by the positional standard deviation o~r, the magnetic 
configurations of the relaxed solutions become increasingly noncollinear. Results are given for the 
magnetic states of a slightly disturbed (ctr/Rq = 0.05), a moderately disturbed (ctr/Ro = 0.10), 
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Figure 1. Illustrations of the magnetic arrangements in two-dimensional nanostructures. The particle 
positions are scattered around the lattice sites of a square array with positional standard deviations 
(a) ctr/Rq = 0.05, (b) ur/Rq = 0.1, and (c) (Jr/Rq = 0.5, where Rq refers to the average interparticle 
distance. The surface coverage amounts to C = 35 %. The left columns refers to microvortex magnetic 
arrangements. Relaxation of the magnetic moment directions using these arrangements as starting 
configurations yields the metastable solutions shown in the right columns. 
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and a strongly disturbed (g r /Rq = 0.50) particle setup having all the same coverage C = 35 %. 
Both the optimal microvortex states and the relaxed solutions are shown. For gr/R = 0.05 the 
MV arrangement resembles quite closely the true solution, see figure 1(a). However, as can be seen 
in figure 1(b), already a moderate amount of positional disorder destroys the MV state. This is 
physically reasonable since the MV order is tightly connected to the presence of a square-lattice 
symmetry of the particle ensemble. The degree of the MV ordering will be quantified below. For 
large g r or for a random particle setup the resulting magnetic arrangement is dominated by the 
formation of chains and loops of magnetic moments with a correlated 'head-to-tail' alignment of the 
particle magnetizations, see figure 1(c) [22]. This reflects the tendency of the dipole interaction to 
favour a locally demagnetized state with a vanishing or small net magnetization. 

In figure 2 we present examples for the dipole energy distribution of metastable states of slightly 
and moderately disturbed particle setups shown in figures 1(a) and (b). These were obtained by 
considering 20000 random initial configurations {0™ ltia1 } and by assigning the energies of the relaxed 
states to energy intervals [E, E + AE], with AE = 5- 10~ 3 K. Distinctly different energy distributions 
are observed for weak and strong disorder. For weak disorder [e.g., gr/Rq = 0.05, see figure 2(a)] an 
asymmetric energy distribution is found. The lowest energy state is reached very often, in fact about 
50 % of the trials relax to that ground state for gr/Rq = 0.03, and about 30 % for gr/R = 0.05. 
Moreover, the ground state configuration resembles closely the MV state with an almost vanishing net 
magnetization. In addition, numerous metastable states are obtained, which energies are distributed 
over a relatively broad range, and which are reached far less frequently. In other words, for weak 
disturbances the catchment area of the ground state is much larger than the ones of the higher-energy 
states. It is interesting to note that the energy differences between the few states with the lowest 
energies are quite larger than those found for less stable magnetic arrangements. Since these latter 
states often exhibit a finite net magnetization, one expects that an external magnetic field tends to 
stabilize these higher-energy metastable configurations. 

Already for moderate disorder g r / R = 0.10 the character of the energy distribution changes 
strongly. An almost symmetric, Gaussian-like energy distribution is obtained around the average 
dipole energy E^ p , see figure 2(b). The number of metastable states has increased remarkedly. In 
fact, out of the considered 20000 random trials no single state is reached twice after relaxation. This 
is in particular true for the low-energy states which are obtained with a very small frequency [see 
the encircled region in figure 2(b)]. For an even stronger disorder the small peak in the frequency 
distribution for the low-energy states disappears completely. Since the total number of metastable 
states increases strongly with increasing disorder, the corresponding catchment areas decreases 
dramatically. The obtained large number of metastable states is consistent with experiments on 
ferrofluids, which show that after application and removal of an external magnetic field the same 
magnetic arrangement is seldomly reached for a second time [23]. 

Furthermore, from figure 2 one observes clearly that the average magnetic dipole energy E dip 
decreases with increasing positional standard deviation gr. In figure 3 we present results for 
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Figure 2. Energy distribution of local energy minima for two disturbed square particle arrays using 
different positional standard deviations (a) <jr/Ro — 0.05, and (b) <jr/Rv = 0.10. The particle 
coverage amounts to C = 35 %. The distributions of the metastable states are obtained from 20000 
randomly chosen initial configurations and sampled into energy intervals with width AE = 5- 10~ 3 K. 
Also indicated is the resulting average dipole energy i?dip(cfl). 



-^dip(°\R) C) per particle as function of o - ^ and for different particle coverages C. A square lattice 
corresponds to a R = 0. Results for random particle setups are also shown. E dip (aii) is calculated 
for the relaxed solutions as well as for the ferromagnetic and microvortex magnetic states averaged 
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Figure 3. (a) Average dipole energy Edi p (<Tn) per particle of a planar array of magnetic particles 
as function of the positional standard deviation tr^. Rq refers to the average interparticle distance, 
and the coverage amounts to C — 35 %. Displayed are Edi p (<JR) for the relaxed magnetic states, 
averaged over 40 different initial conditions (full lines, circles), and the minimum energies (dashed 
lines, squares) which indicate the dispersion of the data. Results are also given for the ferromagnetic 
and microvortex magnetic arrangements, averaged over the in-plane angles. The dotted lines denote 
Edip for a random particle setup with same particle sizes and coverages. All results correspond to the 
average over 20 different realizations of the unit cell. 

(b) Edi P ((JR,C) of relaxed magnetic arrangements as function of <tr for different coverages C as 
indicated. The dotted lines denote Edi P (C) for random particle setups. 

(c) Edi p {<Tn,, C) of relaxed magnetic arrangements as function of the particle coverage C for different 
ctr as indicated. We show -Edip(C) as function of C~ 3/<2 which should yield a linear behaviour from 
a simple scaling estimate. 



over the corresponding in-plane angles 0f m and mv . The obtained minimum values of -Edip(°"i?) are 
also displayed in order to illustrate the spread of the results around the average. Whereas the energy 
distributions shown in figure 2 are determined from a single realization of the unit cell, here E^ p (aR) 
is in addition averaged over 20 different realizations. 
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In figure 3(a), E dip (aji) is shown for a coverage C — 35 %. With increasing o R the average dipole 
energy decreases, i.e., the respective magnetic binding energy increases with increasing disorder [24]. 
For o~r/Rq > 0.5, -Edip(oR) approaches a constant value which corresponds, within the numerical 
dispersion of the data, to the average dipole energy for a random particle array. The decrease of 
the average energy is caused by the nonlinear dependence of the dipole interaction with respect 
to the interparticle distance [24]. In fact, once disorder is introduced, the increase of E dip for 
enlarged distances r^ between some particle pairs is more than counterbalanced by a corresponding 
decrease for smaller distances between other pairs of particles. A similar behaviour is obtained for 
the ferromagnetic arrangement, albeit with a larger E dip (aR). In contrast, for the MV state the 
average energy exhibits a minimum as function of <jr at o~r/Ro ~ 0.15, and approaches -Edip = 
with increasing or. Relaxation in a disordered particle arrangement, see figure 1, is thus crucial to 
the disorder induced reduction of E dip (aR). Let us recall that the FM and MV arrangements usually 
do not correspond to local energy minima for o~r > 0. 

In figure 3(b) we show E d i P (an) of the relaxed solutions and for different particle coverages C. 
First of all one observes that the overall dependence of E dip on or is not significantly affected by 
C. Increasing the interparticle spacing R decreases the magnitude of the average dipole energy, 
which should scale in principle as E dip oc Rq 3 oc C 3//2 . The dependence of E dip on C is depicted 
in figure 3(c) for the relaxed solutions. Indeed, the expected behaviour E d \ p oc C 3 ^ 2 is obtained 
for weak positional disorder <tr. However, for strong disorder and for small coverages a different 
concentration dependence is observed. It seems that the strong magnetic noncollinearities of the 
magnetic configurations render the simple scaling expectation no longer applicable. It would be 
therefore interesting to investigate in detail the energy scaling as function of coverage C especially 
in the limit of strongly disturbed arrays of magnetic particles. 

Note that the magnitude of the dipole energy is comparably small. This is a consequence of the 
disk-shaped particles assumed in our calculations. For compact sphere-shaped particles with the same 
radius r as the disk-shaped ones, thus yielding the same surface coverage C, the corresponding dipole 
energy will be significantly enhanced due to the larger number of atoms N = (ro/ao) 3 per spherical 
particle. Moreover, similar as the positional disorder a particle-size dispersion yields noncollinear 
magnetic arrangements even if the particle centers form a periodic lattice. However, one observes 
that -Edip does not vary strongly with increasing size dispersion [25]. 

Furthermore, we discuss the crossover from a quasi-periodic to a random particle ensemble. In 
figure 4 we show the coverage dependence of the positional standard deviation a R , above which 
the average energy E d \ p {o~R) of a quasi-periodic particle arrangement reaches the limiting value of 
a random setup. One observes that Wr decreases with increasing coverage C, as can qualitatively 
explained by the following simple scaling consideration. The particles are scattered within a certain 
distance around the sites of the square lattice, the average scattering radius R aR can be estimated 
by R aR ~ ctrRq. For a given coverage C the average dipole energy for a random particle setup 
is approximately reached for the crossover standard deviation <7r for which R aR plus the particle 
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of a random particle setup. The symbols refer to the inflection points of -Edi P (o\R) extracted from 
figure 3(b), and the full line to equation 5 using 0.72 as proportionality factor. 



radius equals half of the average interparticle distance, 

Rs R + r ~ R /2 . 
Using C = ir(r /R ) 2 , one obtains 



(4) 



i 

2 



(5) 



In figure 4, <tr(C) is compared with the inflection points of -Edi P (o"i?) extracted from figure 3(b) 
for various coverages C . One observes that except for the largest considered coverage C = 35 % a 
satisfactory agreement is obtained, which supports the validity of the previous scaling considerations. 



Results for the global and local microvortex order parameters M|J° bal (crij) and M^^or) are 
presented in figure 5, assuming two different types of disorder. In figure 5(a) we consider a disturbed 
square particle lattice as function of the positional standard deviation a R for a coverage C = 35 %, 
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Figure 5. Global and local microvortex order parameters M|J° bal {a r) and M^ al (&R) a s function of 
the distortion of the particle array assuming (a) a disturbed square lattice with positional standard 
deviation or and particle coverage C — 35 %, and (b) a square array with randomly distributed 
vacancies with concentration C vac . These order parameters are averaged over many different initial 
arrangements and different realizations of the unit cell, in addition M^ &1 ((jr) is averaged over all 
four-particle plaquettes in the unit cell. The resulting standard deviation is indicated by the error 
bars. 
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and in figure 5(b) we vary the concentration of vacancies in an otherwise periodic array. The rather 
large dispersions of the MV order parameters result from many different initial arrangements and 
different realizations of the unit cell. Locally the MV order is preserved even for strong disorder. 
For example, for o~r/Rq = 0.5 which refers to an almost random particle array, see figure 3(a), we 
obtain M^ al ~ 0.5, and M^ al ~ 0.7 for C vac = 20 %. In contrast, positional disorder has a much 
stronger effect on M|J° bal . Especially positional disturbances quickly destroy the long-range MV 
ordering in the square lattice. For instance, for a particle coverage C = 35 % this occurs already for 
o-r/Ro > 0.10, which is consistent with the minimum of E dip (a R ) of the MV state, see figure 3(a). 

Finally, we present results for the angular distribution of the microvortex angle mv for different 
degrees of disorder. We consider first the effect of vacancies in an otherwise periodic square particle 
lattice, as has been done previously by Prakash and Henley [10]. According to these authors the 
preferred MV angles for a dipole-coupled square spin lattice with a small amount of vacancies 
correspond to the diagonal directions. Indeed, as can be seen from figure 6, a strong frequency 
of the diagonal MV angles mv = 45° and 135° is obtained for low vacancy concentrations C vac . The 
frequencies of these states decrease with increasing C vac , and vanish for C vac > 6 - 7 %. In addition 
to the diagonal states we also obtain a strong frequency of the MV angle along the axial directions 
(columnar states) with MV angles mv = 0°, 90°, and 180°. These states are also present for larger 
vacancy concentrations. We have investigated whether the angular distributions are affected by the 
finite size of the unit cell, in particular whether this could induce the pronounced frequencies of the 
columnar states. If this is true, the ratio of the frequencies of mv = 0° and mv = 90° is affected 
by a variation of the aspect ratio of the rectangular unit cell. We found that for not too small C vac 
the dependence of the calculated frequencies for columnar and diagonal states as function of size and 
shape of the unit cell is quite weak. Hence, our results indicate that the presence of a few vacancies 
induces a more complicated angular distribution of the MV angle. In any case, we would like to point 
out that at present we cannot rule out completely the possibility that the pronounced frequencies of 
the columnar states might be caused by the finite size of the unit cell. 

Concerning other types of disorder, we observe that for a disturbed square lattice characterized 
by the positional standard deviation o~r there is no particular preference of the diagonal states. For 
this type of disorder only the columnar states exhibit a pronounced frequency as compared to the 
rest of the MV angles. The frequency peaks for columnar states decrease for increasing disorder and 
vanish above o-r/Rq ~ 0.20 for the particle coverage C = 35 %. A marked dependence of the angular 
distribution on the size and shape of the unit cell is not observed also for this type of disorder. 
To conclude this discussion, we would like to emphasize that the angular distributions presented 
in figure 6 were all obtained for a nonrelaxed MV state. We have also performed corresponding 
calculations for the relaxed solutions, and have determined the MV states which are closest to 
the noncollinear relaxed magnetic arrangements. In contrast to the results shown in figure 6 the 
corresponding angular histograms for the relaxed states exhibit a strong dependence on the size and 
shape of the unit cell. Therefore, further investigations would be necessary in order to clarify this 
matter. 
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Figure 6. Angular distribution of the microvortex angle </> mv for different vacancy concentrations C vac 
on a square lattice, (a) C vac = 2.0 % (5 vacancies per unit cell), (b) C vac = 3.9 % (10 vacancies), and 
(c) C vac = 5.9 % (15 vacancies). The size of the unit cell is 16 x 16 particles, 1000 different realizations 
of the unit cell are performed. Due to the mirror symmetry we show the angular distribution for 
-20° < mv < 200°. 
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4. Conclusion 

In the present study we have investigated the low-energy properties of disordered planar arrays of 
magnetic nanoparticles interacting by the dipole coupling. Different kinds of disorder have been 
considered. Already small deviations from a square particle lattice lift the continuous degeneracy of 
the microvortex ground state. A strongly noncollinear magnetic order appears in the nanostructure 
which destroys the MV arrangement, as can be seen from the corresponding order parameter. We 
have shown that with an increasing disorder the energy distribution of the metastable states changes 
qualitatively. Our results indicate that it should be very difficult to identify the ground state energy 
and its magnetic configuration for strongly inhomogeneous particle arrangements. This behaviour is 
typical for spin-glass systems or random magnets. A detailed investigation of the spin-glass hehaviour 
of the nanostructure is certainly worthlike. This would require the consideration of additional 
properties such as the nonlinear susceptibility and the time dependent correlation functions [3, 17]. 

It has been shown that the average magnetic dipole energy of an ensemble of magnetic particles 
decreases with increasing positional disorder. Moreover, the average dipole energy Ea\ p of a quasi- 
periodic particle arrangement resembles the one of a random particle setup above the coverage- 
dependent crossover standard deviation a R (C). A simple scaling behaviour of <7r(C) has been 
derived, which reproduces quite accurately the coverage dependence of (7r(C). Structural disorder 
and strong magnetic noncollinearity effects result in a deviation from the straightforward scaling 
i^dip oc C 3 / 2 of the average dipole energy as function of the coverage. 

For a square particle lattice with a small vacancy concentration we found a preferred orientation 
of the microvortex angle mv along the diagonal directions. This is in agreement with the results 
obtained previously by Prakash and Henley [10]. The pronounced frequency of this state decreases 
with an increasing vacancy concentration. However, we also find indications that the columnar states 
are particularly frequent. The difference with previous results could be related to the fact that isolated 
vacancies were considered in [10], whereas in the present study the vacancy concentration cannot be 
chosen arbitrarily small. In contrast to the case of vacancy-induced disorder, for a disturbed square 
particle array characterized by the positional standard deviation <xr a pronounced frequency is found 
only for the columnar but not for the diagonal states. 

The determination of E^ip allows to distinguish whether a magnetic particle ensemble can be 
considered as a weakly or a strongly interacting system. In the first case the interactions can be 
treated as perturbations to the single-particle couplings, whereas in the latter case they have to 
be considered explicitly. The present study can be easily extended in order to take into account 
single-particle anisotropies with distributions of their magnitudes and easy axes. The importance of 
collectively ordered magnetic states for strongly interacting particle systems has been pointed out. 
An increasing average binding energy l-E^pl favours magnetic ordering and should cause an increase 
of its critical spin-glass temperature. In this context the effect of disorder in the particle ensemble 
on the critical temperature is of considerable interest. Also the magnetic relaxation will depend 



Low-energy properties of two-dimensional magnetic nano structures 



17 



sensitively on the degree of disorder in such nanostructured systems. Finite temperature effects 
can be introduced in the framework of a mean field approximation or by performing Monte Carlo 
simulations. Magnetic hysteresis loops and susceptibilities can be derived by applying magnetic fields 
with different directions and strengths. 
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Appendix A. Dipole interaction between magnetic particles with a finite size 

In this Appendix we determine the leading correction to the magnetic dipole interaction, see 
equation 1, beyond the point-dipole sum. The finite extension of the particles for the most general 
case of arbitrary sizes and shapes is taken into account. For a hexagonal lattice the finite particle 
size has been considered already by Politi and Pini [12]. First, the interaction between a particle 
pair is calculated. Then a particle ensemble with an infinite lateral extension is modeled by means 
of a rectangular unit cell with periodic boundary conditions. 

The dipole- quadrupole correction 

Consider a particle % containing Aj atoms on lattice sites k with position vectors = + r k . The 
center of gravity of this particle is given by 



The finite size of particle i is taken into account by the quadratic deviations with respect to its 
center, 



and similarly for (yf) and (zf). These quantities depend on the size and shape of the particle and 
have the dimension of an area. Furthermore, we define by 




(A.l) 





r ik = ( x i0 + x k f + (ViO + Vkf + (z i0 + z k f 



rj (l + e ik ), 



(A.3) 
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e ik = \ ( 2 x io + 2 Vio Vk + 2 z i0 z k + x 2 k + y 2 k + zf\ , 
r iO v 7 



(A.4) 



with rf = x 2 + y 2 + z 2 . The factor (1 + e ik ) appearing in the denominator of the dipole interaction 
energy, see equation 1, is expanded to second order in x k , y k , and z k as 



1 



. 5 35 



(A.5) 



(l + e tfe ) 5 / 2 

1 - ^-2- (2 x i0 Xfc + 2 y i0 J/fc + 2 ^ + ^ + ^ + ^) + — — (x i0 x k + y i0 y k + 2*) • 

Now all different sums over the atomic sites k of particle % are performed up to this order. The 
non- vanishing terms are listed in the following equations A.6-A.15. 
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(A.13) 



(A. 14) 



(A.15) 



A simple extension of the previous considerations allows to calculate to the same order the 
corresponding sums involved in the interaction between two extended particles % and j with sizes 
Aj and Nj. Let r^o = (a^o, 2/ijo? %o) denote the relative position vector between the particle 
centers. The non- vanishing sums in the interaction energy are given by equations A.6-A.15 after the 
replacements Nj — > NNj, rjo — > r^o, xjo — > x^o, (a^) - ► (a^ 2 ) + (a^ 2 ), etc., have been performed. 
Concerning the summation over all atoms k' of particle j, one obtains, for example, for equation A. 6, 
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(A.16) 



Thus, within this expansion the dipole interaction between particles % and j is expressed in terms of 
the distance r^o between their centers and the quadratic deviations (xf), (a; 2 ), etc., characterizing 
the sizes and shapes of the two particles. The point-dipole sum is recovered by setting (xf) = 0, 
etc. The correction to the point-dipole sum is a dipole-quadrupole interaction, being of the order 
((xf) + {Xj))/rfj , i.e., the ratio of the particle extensions and the square of the interparticle distance. 
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distance d / a Q 

Figure Al. Comparison between the point-dipole approximation, the point-dipole approximation 
augmented by the dipole-quadrupole correction, and the exactly calculated dipole interaction E dlp (d). 
We consider two square-shaped flat (001) magnetic particles, consisting of 11 x 11 = 121 atomic 
magnetic moments each. Ea p (d), in units of E® ip = /U^/ag, is shown as function of the distance 
d between the particle centers in units of the lattice constant a , where /z a t is the atomic magnetic 
moment. Note that for d/ao = 10 the particles touch each other. The thick arrows illustrate the 
magnetization directions, which are oriented head-to-tail. 

Evidently, the effect of finite particle sizes becomes more important the closer the particles are located. 
This is the case, for example, for a densely packed 3D ferrofluid or for a layered nanostructured 
particle ensemble with a large surface coverage. 

In figure 7 we compare the exact dipole interaction energy E^ ip with the results obtained by 
using the point-dipole approximation with and without the dipole-quadrupole correction. E^ p (d) 
for two square-shaped particles is shown as function of the distance d between their centers. Notice 
that the correction for nearby particles amounts to ~ 25 % of the point-dipole sum in this case. For 
simple particle geometries (e.g., disks, cubes, or spheres) the quadratic deviations (xf), etc., can be 
calculated analytically. In case of spheres the dipole-quadrupole correction vanishes since a sphere 
has no quadrupole moment [26]. 
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Let us now consider a planar particle ensemble infinitely extended in the xy plane. The unit cell 
has the size L x x L y and consists n particles. Periodic boundary conditions are applied laterally, 
whereas the vertical extension along the z-direction is finite. The dipole interaction sum runs over all 
particle pairs within the same and between different unit cells. The position vectors connecting two 
particle centers are given by r^o = (xijo + l x L x , y^ + l y L y , %o), with l x and l y integers. For such a 
periodic planar system an Ewald summation technique can be applied by using a rapidly convergent 
2D lattice summation [20] . The following general sums over the unit cells need to be considered, 



rpa/3-/ 



VijO 



(A.17) 



ijO 



The prime at the sum indicates that the term with r^o — is omitted. The different lattice sums 
appearing in equations A.6-A.15 can be obtained by an appropriate choice of the positive integers 
a, (3, and 7. In general one has to distinguish between the cases 2^0 = and z^q 7^ 0. 

The dipole field Bf ip of particle i, see equation 1, has the components 
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(sin^j cos 4>j , sin 9j sin (f>j, cos 0j) determines the 



direction of the particle magnetic moment, with 9j and <f>j the polar and azimuthal angles. 

Introducing; now the lattice sums T t f 7 and after some algebra one obtains for the three 
components of the dipole field 
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(A.21) 



(A.22) 



(A.23) 



The lateral extensions L x and L y of the unit cell and the components and y^Q of the 
interparticle distances within a unit cell are involved in a complicated manner in the lattice sums , 
whereas the dependence on the vertical distances z^ appears explicitly. Note that the contributions 
from mirror particles i = j located in different unit cells are taken into account. The expressions 
for the dipole field components, equations A.21 - A.23, are simplified appreciably if all particles are 
located in the same plane (%o = 0). In this case B*' dip oc m|, hence the z component of the dipole 
field vanishes for an in-plane magnetization (m| = 0). As before, the point-dipole sum is recovered by 
setting (xf) = 0, etc. The dipole energy E^ip per unit cell is obtained from equation 1 by performing 
the sum over all particles i — 1 . . . n in the unit cell. 
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